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Abstract 

We compute the free energy density F for gauge theories, with fermions, at high tem¬ 
perature and zero chemical potential. In the expansion F = T 4 [co + C 2 g 2 + c%g 3 + (<f 4 In g + 
cf)g A + (c' 5 In g + c 5 )g 5 + 0(g 6 )], we determine c' 5 and c 5 analytically by calculating two- and 
three-loop diagrams. The g 5 term constitutes the first correction to the g 3 term and is for 
the non-Abelian case the last power of g that can be computed within perturbation theory. 
We find that the g 5 term receives no contributions from overlapping double-frequency sums 
and that c' b vanishes. 





I. INTRODUCTION 


The perturbative expansion of the free energy density of high-temperature gauge theory 
in four dimensions can be written as 

F = T 4 [co + C 2 <T + c 3 g 3 + (c' A \i\g + c±)g l + (c' 5 \ng + c§)g b + 0(g 6 )}, (1.1) 

where the care numerical coefficients (that depend on the field content of the theory, 
the renormalization scheme and the renormalization scale) and where we have assumed 
the temperature high enough that fermion masses can be ignored. Previously, F has been 
computed to 0(g 3 ) for QED by Akhiezer and Peletminskii [|l| and for QCD by Kapusta ||, 
while the g 4 In g term was obtained by Toimela || . More recently, F has been computed 
to 0(g A ) by Arnold and one of the present authors (C.Z.) Q. Coriano and Parwani || 
have recently studied high-temperature QED up to 0(g 5 ). The free energy density (or, 
equivalently, the pressure) is also known to 0(g 5 ) in $ 4 theory (see Ref. || and references 
therein). Here we determine the coefficients c ' 5 and C 5 in expansion (|1 . 1| ) . For this purpose we 
need to take into account Debye screening at three loops (for a review on Debye screening, 
see Refs. [7]||). 

Note that, for the non-Abelian case, the g 5 term is believed to be the last power in g 
accessible within perturbation theory || (for a review, see Refs. fj],[J). Starting at four loops, 
infrared problems that are believed to be cured by non-perturbative magnetic screening lead 
to contributions to the g 6 term from diagrams with arbitrarily high numbers of loops. 

The g 5 term is also interesting because it constitutes the first correction to the g 3 term, 
the lowest order at which Debye screening plays a role. The renormalization group invariance 
to this order can be tested. The dependence of F on the renormalization scale due to the g 3 
term should be diminished by including the g 5 term. Checking this, we can gain some idea 
about the theoretical uncertainties of the g 3 term as well as the behavior of the perturbative 
expansion. Also, our result can be used for a test of an evaluation of F on the lattice. 
Finally, our result is potentially interesting for the evolution of the early Universe, where 
one might have to add scalars to the theory. 
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In section [TJ notation and conventions are established. In section 0 we outline our gen¬ 
eral computational procedure and emphasize what is new as compared to the g A calculation. 
In section we conclude by presenting and analyzing our result as well as comparing it to 
related results. 


II. NOTATION AND CONVENTIONS 

We use the same notation and conventions as in Ref. J4|. We now present an almost 
verbatim review of these to keep this work as self-contained as possible. 

We consider gauge theories given in Euclidean spacetime by Lagrangians of the form 

Ce = H (% - 4 , + i (d„Al - d„A“ + gf^A^Atf + £ gt , gh . (2,1) 

with gauge fixing and ghost term £ g f jgh , and where the T a are the generators of a single, 
simple Lie group, such as U(l) or SU(3). To simplify our presentation, we will not derive 
results for an arbitrary product of simple Lie groups such as SU(2)xU(l), but such cases 
could easily be handled by adjusting the overall group and coupling factors on the results 
we give for individual diagrams. 

d A and Ca are the dimension and quadratic Casimir invariant of the adjoint representa¬ 
tion, with 

S aa = d A , fabcjdbc = CAS ad (2.2) 

dp is the dimension of the total fermion representation ( e.g ., 18 for six-flavor QCD), and S'f 
and S* 2 F are defined in terms of the generators T a for the total fermion representation as 

S F = -y-tr(T 2 ), S 2 f = J-tr[(T 2 ) 2 ], (2.3) 

«A «A 

where T 2 = T a T a . For SU(1V) with ri{ fermions in the fundamental representation, the 
standard normalization of the coupling gives 

1 TV 2 —1 

d A = N — 1, Ca = TV, dp = TVnf, S'f = -rif , S 2 p = ^ ri{. (2.4) 
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For U(l) theory, relabel g as e and let the charges of the m fermions be q^e. Then 


d A = l, C A = 0, d F = n { , S F = J2Vi, S 2F = J2 ( ii- (2-5) 

i i 

If the fermion representation is irreducible or consists of several identical copies of an irre¬ 
ducible representation [as in Q2.4| ) above], we have 

d A S F = d F S 2 F . (2.6) 


We work in Feynman gauge. We also work exclusively in the Euclidean (imaginary time) 
formulation of thermal held theory. We conventionally refer to four-momenta with capital 
letters K and to their components with lower-case letters: K = (fc 0 , k). All four-momenta are 
Euclidean with discrete frequencies ko = 2nnT for bosons and ghosts and ko = 2 tt T 

for fermions. We regularize the theory by working in d = 4—2e dimensions with the modified 
minimal subtraction (MS) scheme, which corresponds to doing minimal subtraction (MS) 
and then changing the MS scale /i to the MS scale jl by the substitution 


h 2 = 


e 7E /x 2 

47T 


(2.7) 


The trace over the identity in spinor space is by convention tr/ = 4. 

To denote summation over discrete loop frequencies and integration over loop three- 
momenta, we use the shorthand notation 


for bosonic momenta and 



d 3 2e p 

(27r) 3_2e 



for fermionic momenta, where 


( 2 . 8 ) 


(2.9) 


Y. — Y. . Y. — Y ■ (2.io) 

PO p 0 =2irnT {p 0 } p 0 =2n (n+±)T 
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FIG. 1. The one-loop gluon self-energy. 


We handle the resummation of hard thermal loops [which is required to make pertur¬ 
bation theory well behaved beyond 0(g 2 )\ as was done in Ref. Q. Specifically, we must 
improve our propagators by incorporating the Debye screening mass M for which is 
determined at leading order by the self-energy diagrams of Fig. |1]: 


M 2 S ab = ngg(O) = n"( 0 ) = g‘8 


jab 


^2?ab 


c A 


{d - 2 ) t^- is ^ d - 2 >Lw 


1 

Xq}Q 2 


( 2 . 11 ) 


This is accomplished by rewriting our Lagrangian density, in frequency space, as 


C E = (£ e + i«’«(,,) - ±M 2 A‘A‘S Po , 


( 2 . 12 ) 


where 5 Po is shorthand for the Kronecker delta symbol h Po ,o- Then we absorb the first A q 
term into our unperturbed Lagrangian Cq and treat the second A q term as a perturbation. 

Since the free energy density is computed by considering vacuum diagrams (diagrams 
without external legs), there is no need to explicitly introduce wave function renormaliza¬ 
tions. We only need to renormalize the coupling constant. We do this by expressing the 
bare coupling constant g b in terms of the renormalized coupling g, 

/ 11C a -45 f g 2 x 


= z y = 1 


+ 0(pW, 


(2.13) 


v 3e (47t) 2 j 

and then using the bare coupling constant g b to the required order in g at all vertices in the 
vacuum diagrams. Through g 5 it is sufficient to know the one-loop renormalization given 
above. However, for the computation of the Debye screening mass fl2.11| ) we have used 
the renormalized coupling g. This is allowed because for the cure of the infrared problems 
achieved by reorganizing the perturbation expansion according to (|2.12|) , only the leading 
contribution in g to the Debye mass is crucial. 
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III. COMPUTATIONAL PROCEDURE 


A. g 5 Order Contributions 

In the expansion of the free energy density, the zeroth-order term represents the free 
energy density of an ideal gas containing free gauge bosons and massless quarks. The 
leading contribution due to the interaction is of order g 2 which is represented by two-loop 
diagrams. For the calculation of higher-order terms the resummation (|2.12|) for the static 
timelike gluon propagator is required because of infrared divergences. It is this resummation, 
which introduces the Debye screening mass into the theory, that causes the expansion of the 
free energy density being in powers of g instead of g 2 . Consequently, one cannot determine 
the order of a diagram by naively counting the number of interaction vertices. The leading 
odd-power contribution is of order g 3 which comes from the one-loop diagram with the 
resummed static gluon propagator. The g 4 term receives contributions from the subleading 
pieces of two-loop diagrams as well as the leading pieces of three-loop diagrams. To get 
the g 5 term of the free energy density, we need to compute the two-loop diagrams to higher 
order and the subleading pieces of three-loop diagrams. Fig. contains all the diagrams 
contributing to the free energy density up to g 5 order. 

We mainly follow the formal manipulations in Refs. J3J to simplify the sum-integrals 
obtained by applying the Feynman rules to the diagrams in Fig. j|. We shall focus on the 
contributions at order g 5 . Due to the resummation, there are two momentum scales, gT 
(the Debye mass) and T, appearing in the sum-integrals. We will conveniently refer to 
momenta of order gT as “soft” and momenta of order T as “hard.” A sum-integral usually 
gains its dominant piece from a momentum integral region where some momenta are hard 
and others are soft. Our strategy is to identify the soft and hard momenta in the sum- 
integrals and carry out an expansion about the ratio between the Debye mass or the soft 
momenta and the hard momenta to extract the leading piece. We then construct a new 
sum-integral from the original one by subtracting its leading piece out. We then find the 
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= - n, v + 'VN^^W 

FIG. 2. Diagrams contributing to the free energy of gauge theories with fermions. The crosses 


are the “thermal counterterms” arising from the last term of ( 2.12 ) 
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corresponding soft and hard momenta for this new sum-integral to extract its leading piece 
which is the next-to-leading piece for the original integral. This procedure enables us to get a 
systematic expansion in powers of g. Since when getting the leading piece of a sum-integral, 
a soft momentum is always neglected compared to a hard momentum, the integrals for two 
different scales T and gT are separated or decoupled. 

In the following two subsections, we give more details about our computational procedure 
by considering separately two- and three-loop diagrams. We concentrate on explaining the 
way in which we perform the calculation but spare the reader from all the messy details for 
computing individual diagrams. The expressions for all the diagrams in Fig. are provided 
in Appendix |Aj. 


B. Two-loop Diagrams 


To illustrate the general discussion in the previous subsection, let us first consider a 
typical two-loop sum-integral arising from the setting sun diagram (e) in Fig. |2|. 

Atyp = fp Q P 2 + M 2 Q 2 (P + Q) 2 ' ^ 

We shall use superscripts to denote the pieces in the expansion of A typ in M, i.e., A{; y p for 
the leading piece, A[ yp for the subleading piece and so on. The leading piece of A[ yp may be 
obtained by setting M to zero since A typ is infrared safe as M —> 0 and since the typical P 
contributing to A typ is of order T which is much larger than M. We can then subtract this 
leading piece from A typ which gives 


A _ /t(°) _ l\/[2 
^typ ~ Vhyp M 


J P 0 


<?o 2 


PQ P 2 (P 2 + M 2 ) Q 2 (P + QY 


(3.2) 


with 


AO) _ -rf fipo % 
typ Tpq P 2 Q 2 (P + Q) 2 ' 


(3.3) 


Now, we find that the P integral in the second term on the right-hand side of Eq. is 
infrared sensitive to M which means this P integral picks up contribution mainly at the 
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region where P is of order M. Thus, P is soft for the subleading piece of A typ . Since Q is 


always hard, the subleading piece A^, is 


d (1) = 

^typ 


-M 2 ^t --- 

Tpq P 2 (P 2 + M 2 ) Q 4 


% 2 = T 4 A 

la Tq g 4 


(3.4) 


and then 


A 


typ 


/t(°) _i_ /IC 1 ) i M 2s >f _ — _ 

typ + typ + Tpq P 2 (P 2 + M 2 ) 


% 2 


% 2 


g 4 g 2 (p + g) 2 


(3.5) 


Here we have defined the integral 


Jla f P 6po {p 2 + M 2 P 2 ) M t P P 2 (P 2?, + M 2 ) ^ 3 ‘ 6 ^ 

which is of order g in four spacetime dimensions. Its value can be found in Appendix |B]. 

Now since the P integral in the last term of ( |3.5| ) behaves as 1/(P 2 + M 2 ) when P is 
much less than T, we find that the P integral receives its main contribution from the region 
where P is hard. Thus, 


and 


A (2) = 


M 2 


xf ^P°_ 

Tpq P 4 



g 2 (p + g) 2 


(3.7) 


|( 0 ) 


/t = A { ) A y } 4 - A [ J 

typ 41f vr) r 41-j- iVr) i 4±f vn 1V1 


x typ 


(i) 

typ 


i (2) 


J P o 


1 typ 


fpq P 4 (P 2 +M 2 ) 


g 4 g 2 (p+g) 2 


(3.8) 


We then identify P as soft to get as 

=4 m2j 4 q |- < 3 - 9 > 

where we have expanded the denominator (P + Q) 2 in powers of the ratio of P and Q, 
replaced S Po (P ■ Q) 2 by S Po p 2 q 2 /(d — 1) and integrated by parts in q. It is not hard to see 
that At°p, A^p, A[ yp, and At yp are of orders 1, g, g 2 , and g 3 , respectively, in four spacetime 
dimensions. 

For a two-loop diagram, there is an extra factor g\ = p~ 2e Z 2 g 2 multiplying A typ . Thus, 
A^p and A^ p contribute to the g 5 part of F. Therefore, the above asymptotic expansion 
of A typ demonstrates how we extract the g 5 order contributions to the free energy density 
from two-loop diagrams. 


d (3) 


7\/f4 


i r-L 


p 2 4 (p • g) 2 

ge 
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C. Three-loop Diagrams 


We now turn to the three-loop diagrams. Previously ||], the Debye mass in the resummed 
propagator (2 . 12| ) has been ignored for the three-loop diagrams since after reorganizing 
perturbation theory the three-loop diagrams are infrared finite if we set the Debye mass 
to zero. Now, since we need to explore one more order, the subleading terms need to be 
extracted. Since the Debye mass only appears in the static gluon propagator and is only 
probed by soft momenta of order gT , the sum-integrals that we need to deal with contain 
at most two frequency sums and at least one soft-momentum integral. 

Three different cases appear: 

(1) Case one is where we have pure three-dimensional triple-momentum integrals with 
the Debye mass being the only mass scale. Since there are three loops, there is a prefactor 
g 4 T' i for these triple-momentum integrals. Thus, these three-dimensional triple-momentum 
integrals will give a result proportional to the Debye mass to make up for the missing mass 
dimension, i.e., they contribute to the free energy density at order g 4 T 3 M oc g 5 T 4 . These 
pure three-dimensional three-loop integrals are one of the new features which we encounter 
in the (f order calculation. Their evaluations are provided in Appendix [C]. 

(2) The second case we consider is where the sum-integral contains only one sum and 
two three-dimensional soft-momentum integrals. Thus, in this case there is only one hard 
loop momentum integral. Neglecting the soft momentum relative to the hard momentum 
enables us to decouple the three-loop integral into a product of a one-loop sum-integral and 
a two-loop pure three-dimensional momentum integral, which can be evaluated by standard 
methods. In fact, we can show that these do not contribute to the free energy density 
at order g 5 . The point is that a two-loop three-dimensional momentum integral produces a 
result proportional to the Debye mass to an even power. Since only soft-momentum integrals 
generate odd powers in g, this case is not relevant to the g 5 order evaluation. 

(3) Finally, we need to consider the case where two sum-integrals and one soft-momentum 
integral are involved. Again, for the hard-momentum sum-integral, it is valid to neglect the 
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soft momentum compared to the hard one. This leads to a product of a two-loop sum-integral 
and a one-loop three-dimensional momentum integral. With the methods developed in Ref. 
|4j], we can evaluate the two-loop sum-integrals. However, it turns out that after we sum up 
all the pieces contributing to free energy density at g 5 order, all the overlapping double sum- 
integrals cancel and only the non-overlapping double sum-integrals, which can be written as 
a product of two one-loop sum-integrals, survive. This observation was already made earlier 
for the case of QED ||. 

As a concrete example, let us consider the simplest three-loop diagram, the basketball 
diagram (j) in Fig. ^]. Applying the Feynman rules gives 




x 


'1-5, 


<70 


Q 7 


+ 


5, 


<70 


- 1 - 


PQK 

1 


'l—5 


^P0 \ 

P 2 ' P 2 +M 2 J 


+ 


A^ 2 (P+Q+A^) 2 


(3.10) 


Q 2 +M 2 J P 2 Q 2 

The definition of J^ u may be found in Appendix [B|. It is convenient to rewrite the expression 
above as 


rbb 


^F, = -(d-l)g 4 d A CU-dlU 


PQK 


5po5<jo5fc 0 


1 


+2 


J v 0 


(P 2 +M 2 )(Q 2 +M 2 ) P 2 Q 2 J K 2 (P+Q+K) 2 

l-5fc 0 


(3.11) 


( S P0 


( ^<70 

^<70 \ 

\P 2 +M 2 

P 2 J 

Iq 2 +m 2 

Q 2 


5 P0 \ 


1 - 5„„5 


qo°ko 


fpQK \P 2 +M 2 P 2 J Q 2 K 2 (P+Q+K) 2 J ' 

Let us consider each term at the right-hand side of Eq. (|3.11|) above. The first term, involving 
J^n, is of order g 4 and represents the leading contribution of many three-loop diagrams and 
has been evaluated in Ref. Q . 

Recall that we listed three cases in the general treatment of the three-loop diagrams. 
The second, third, and fourth terms correspond to these three cases, respectively: 

The second term is a three-loop pure three-dimensional momentum integral. This is the 
first case discussed above. We encounter a class of these three-dimensional integrals which 
are defined in Appendix [TJ and computed in Appendix 0. 
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The third term involves two soft-momentum integrals corresponding to the second case. 


Since K is hard, it is valid to neglect P, Q compared to K to write this third term as 


S K S K 

o J PQ 


P' 2 +M 2 


1 

P 2 


Q 2 +M 2 Q 2 I Tk K a 


3(d — 1) 




(3.12) 


which is of order g e as what we have expected (even power in g ). 

The fourth term corresponds to the third case where P is the soft momentum since P 
needs to resolve the Debye mass M. Therefore, as described above, we can approximate 
(P + Q + K) 2 as (Q + K) 2 (the case where q 0 + k 0 = 0 and \q + k\ is soft contributes only at 
order g (j because of phase space suppression). Now, the P integral decouples from the Q , K 
sum-integral. Therefore the fourth term in (|3.11 ) can be expressed as a product of a single 
three-dimensional momentum integral and a double sum-integral: 


1) 4j j jb 


(3.13) 


where we have introduced as 


I b = i - 

sun TqkQ 2 K 2 {Q + K) 2 


(3.14) 


and omitted the term vanishing in dimensional regularization. 

Therefore, we have explicitly shown how to extract the order g 5 contribution to the free 
energy density from a simple three-loop diagram. These order g 5 contributions are expressed 
as either three-loop pure three-dimensional momentum integrals or as products of a single 
three-dimensional soft-momentum integral and a double sum-integral. For other three-loop 
diagrams, parallel steps can be followed except for possibly more elaborate expressions; of 
course, we need to introduce more double sum-integrals and pure three-dimensional momen¬ 
tum integrals, ffowever, as mentioned before, when we add up all the pieces contributing to 
the free energy at order g 5 from the three-loop diagrams, these double sum-integrals cancel 
except for the “non-overlapping” sum-integrals that can be expressed as a product of two 
single sum-integrals. We do not know a fundamental reason which leads to this cancellation. 
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IV. RESULT AND ANALYSIS 


Combining the results for all the diagrams listed in Appendix [A] produces the final result 
for the free energy density through order g 5 in four spacetime dimensions as 


F = d A T 4 


71 

~9~ 


1 (, 7 d F \ 

5 ( 1 + id!) + V4W 




-48 


9_ 

A-K 


3 /Ca + ^f\ 3/2 


—y 

47T/ 


_ 48( £)A A (C A + SF ) ln UJ<±±* 


C, U' n i«T + 3 C(—3) 3 C(—1) 7 ' + 5 , 




1759 37, 

et— —I—— hi 2 


F \ 3 4ttT + 3C(-3) 3 C(—1) 7E 3 + 5 


60 5 

In 2 


105 


+S2F 


£) 5 m I/2 


-24 In 2^ 


C\ (176 In-^-+1767 e —247t 2 —494+264 In 2 ) 


4v tT 


+C A S F (112 In yv 


f! 127 E + 72-128 In 2 


+S 7 ( —64 In 




-144 A- 


2F 


4ttT 
+ 0(g*) 


-647e+ 32—128 In 2 


(4.1) 


where ( is Riemann’s zeta function and 7 B is the Euler-Mascheroni constant 
For QCD with nj quark flavors, to g 5 order, the free energy density is 

8 vr 2 T 4 


F = 


45 


+g 4 


1 + §n f - 0.09499 g 2 (l + ^n f ) + 0.12094 g 3 (l + ±n f 
0.08662 (l + ±n f ) In ( gy J 1 + |n f ] - 0.01323 (l + ^n f ) (l - |n f ) In ^ 


+0.01733 - 0.00763 n f - 0.00088 n\ 


+<7v 1 + \ n i 




0.02527 (l + \m) (l - |n f ) In ^ 

-0.12806 - 0.00717U/ + 0.00027rrJ +0(/ s 


(4.2) 


where we have evaluated the coefficients numerically. 
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For QED with rif charged fermions with charges q t e, the fifth-order free energy density 
may be read from the expression above on using (|2.5| ) 

^.2^4 y v 1 /o / o \ 5 


C">(5) _ 

^QED — 


9^ 


(Ei 2 


i\ l,2 (±-X 

Air) 


(E 9?) ' ( 1 84 In + 647 e - 32 + 128 In 2) + 144 £ q, 


(4.3) 


Taking g — T, it is not hard to check that this result agrees with Ref. |SJ]. 

As in Ref. P|, we now check whether the perturbative expansion of the QCD free energy 
density behaves well for physically realized values of couplings to g 5 order. Although the 
free energy does not have any renormalization scale g dependence, the partial sum does so. 
If the perturbative expansion is well behaved, including higher-order corrections into the 
partial sum reduces the /i dependence. The inclusion of the g 5 order term in the partial sum 
compensates the g dependence of the g 3 teruf]. Besides looking at the g dependence of the 
partial sums, we also compare the size of the contributions from each order. 

Define a s (T) = g 2 (T)/(Air). Fig. [| shows the result for six-flavor QCD when aq(T)=0.1 
(which corresponds to scales of order a few 100 GeV). The free energy density is plotted vs 
the choice of renormalization scale g. We have taken 

m * m - A ‘4 + 1 111 (* - ^ T) ln t) ■ < 4 - 4) 

where 

/?0= (4§F + N • A = (ijl (-f Cl + f CaS f + 8S 2F ) . (4.5) 

In Ref. [(J], it was found that including the g 4 term does not make the partial sum for the 
free energy density less dependent on the renormalization scale. There, one of the main 
sources of the g dependence is the g 3 term which requires the order g 5 term to balance its 


*In Ref. Q, with the help of the renormalization group a g r=> ln g/T was introduced to compensate 
for the g dependence due to the g 3 term. We note that our present g 5 order result correctly 
produces this desired g 5 ln g/T term. 
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FIG. 3. The dependence of the free energy density F on the choice of renormalization scale Jl 
for six-flavor QCD with a s (T ) = 0.1. The free energy density is normalized in units of the ideal 
gas result — + j^odp)Tr 2 T A . The dotted, dot-dashed, dashed, and solid lines are the results 

for F including terms through orders g 2 , g 3 , < 7 , and g 5 , respectively. 


renormalization scale dependence. According to Fig. [3], inclusion of the g 5 term in the partial 
sum does not generally make this sum less dependent on /i and the perturbative expansion 
does not behave well in this respect. For Jl — T, the terms at each order are 

79v t 2 T 4 


F = 


90 


1 - 0.0846 + 0.0976 + (0.0255 + 0 - 0.0192) + (0 - 0.0818) + 0{g b ) . (4.6) 


For this value of a s , the g 2 and g 3 terms have about the same size. This does not necessarily 
mean that perturbation theory does not work well since the g 3 term is the leading term of 
new physics at the scale of gT instead of being a correction to the g 2 term. If the corrections 
at g A and g 5 are smaller than the g 2 and g 3 terms, perturbation theory may still work well. 
However, the numerical values above show that the g 5 term appears not to be generally 
smaller than the g 3 term. Therefore, perturbation theory seems not to work well for this 
value of a s which corresponds to QCD at the electroweak scale. 

Fig. f| shows the /i dependence for a s (T) = 0.02 without any fermions, he., with rif = 0. 
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i°g 10 (M-/T) 


FIG. 4. The same as Fig. |3] but for a s (T) = 0.02 without fermions. 


This is interesting since it is [see Eq. (|2.4|) 1 equivalent to pure SU(2), he., electroweak gauge 
theory at the electroweak scale, with ct w ~ 1/33. It is not hard to see that the g 5 order free 
energy density is less sensitive to the renormalization scale than the g 4 order free energy 
density. However, the free energy density through g 5 order is not more stable than the result 
through g 3 order. This is due to a large cancellation between the g 3 term and the g 2 term. 
Here are the values for the contributions at each order to the free energy density for the 
choice jj, = T: 
d A 7T 2 T 4 


F = 


45 


1 - 0.0239 + 0.0152 + (-0.00378 + 0 + 0.00109) + (0 - 0.00406) + 0(g 6 ) . 

(4.7) 


Obviously, the corrections at orders g 4 and g 5 are smaller than the g 2 and g 3 terms. This 
suggests that perturbation theory works. 

In Fig. || we show a similar plot for rif = 6 and a s (T ) = 0.001 where the behavior of the 
perturbative expansion is good. In Fig. ||. we provide the corresponding plot for a s (T ) = 0.2 
but rif = 5 (for T being several GeV). 

We like to comment on the absence of the g° in g term in the expansion of the free energy 
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FIG. 5. The same as Fig. § but for a s (T) = 0.001. 



FIG. 6. The same as Fig. j| but for a s (T ) = 0.2 and rif = 5. 
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density. It is convenient to view the contributions to the free energy at each order with 
an effective held theory technique as was done in Refs. [lOj. In hot gauge theories, there 
are three relevant scales in the imaginary time formalism: T, gT , and g 2 T. The scale T is 
related to the nonstatic fields while gT is the scale for the Debye screening effect. The scale 
g 2 T is believed to be the inverse of the magnetic screening length which cures the remaining 
infrared problem of hot non-Abelian gauge theories. Since the g 2 T scale contributes to 
the free energy starting only at order g 6 , we can ignore it. Imagine first integrating out 
the nonstatic fields (scale T physics) to arrive at an effective held theory which correctly 
describes physics in the low energy region (of order gT). Let A be the cutoff separating the 
scales T and gT. Integrating out the nonstatic helds gives a contribution to the free energy 
density f(T, A, g 2 ) which has an even power expansion in g since no resummation is required. 
Integrating out these nonstatic helds also generates effective interaction terms for the static 
helds. This introduces a dependence on the cutoff A into the bare parameters of the effective 
held theory for the static helds at scale gT with cutoff A. An lng term arises only through 
the logarithm of the ratio of the scales gT and T, i.e., through cancellation between In A /T 
and In A/(gT). In A /T terms enter the free energy density through f(T,A,g 2 ) and the bare 
parameters of the effective theory for the static held at scale gT. It can be shown that there 
are only two parameters relevant to the free energy density through g 5 , the effective mass and 


the effective coupling constant |l(|. Since the coupling constant in the (superrenormalizable) 
effective theory requires no renormalization, there will be no In A /T appearing inside the 
effective coupling. Since /(T, A,g 2 ) has an even power expansion in g, a g 5 lng term comes 
only from the cancellation between an In A /T term in the effective mass term and another 
In A/ (gT) coming from the evaluations of the effective theory. Therefore, the absence of 
the g 5 lng term means that there is no InA/T at order g 4 in the effective mass term which 
has been examined explicitly in Ref. [|TT[ . In other words, this effective mass has vanishing 
anomalous dimension and does not “run” at the leading order as we vary the cutoff A [|T2| . 
In fact, at the next-to-leading order, the effective mass does “run” |T3 . 


As an outlook, it would be interesting to investigate the reasons for the cancellation 
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of overlapping double-frequency sums as well as for the absence of a g 5 In g term in F. It 
would further be worthwhile to include scalar fields and to consider the case of non-vanishing 
chemical potential. 


Note added: Recently, using a different method, Braaten and Nieto have confirmed our 
result (ITS . 


We thank T. Clark and S. Love for useful discussions and advice. We are grateful to 
P. Arnold for suggestions on analyzing the behavior of perturbation theory. We also thank 
E. Braaten and A. Nieto for beneficial discussions and communications. This work was 
supported by the U.S. Department of Energy, contract No. DE-FG02-91ER40681 (Task B). 
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APPENDIX A: RESULTS FOR INDIVIDUAL DIAGRAMS 


Here are the contributions to the free energy density through g 5 order from individual 
diagrams in Fig. ^]. The newly appearing symbols are defined in Appendix [Bj. 


- /i 2e F a = -\d A df In P 2 - d A -l-M 2 J la (Al) 

2 Jp d — 1 

-g 2e F b = d A ^lnP 2 (A2) 

-/i 2e F c = l -d A M 2 J lA (A3) 

-^F d = -?itd A C A (d - 1 ){db\ + 2AA) (A4) 

Z 2 n 2 

~^F e = —^-d A C A {—18(d - 1)6 2 + 6(d - l)(2d - 7 )J la b l 
+6[J 2 a + 4M 2 J 2a + 4 M 2 A\ - 2(2 d - 3)M 2 A b ] 

+(2d 2 -13d + 39)M 2 J la b 2 } + 0(g e ) (A5) 


= -^fd A C A [6b 2 - 6(d - 3)Jigbi + 12M 2 A b - (d - 5)M 2 J la & 2 ] + 0(g 6 ) (A6) 

= f^ 2 {(20d - 23)/ b b a b u 

+4J la [4/ s b un + (16d - 23)/ 2a - (16d - 28)/ 2b - 2(4d - 7)(d - 3)M 2 ] 


-^Fi 

F) 

-^F k 

-^F t 


+2(3J 3a + 4J 3b + 16J 3c — 36J 3e + 16J 3 i + 12 J 3 j — 12J la J 2a )} + 0(g 6 ) 


<V_ 

16 


d A C‘i {-/ bb u + 4J la [—/ 2a + 2 J 2b + (d - 3)M 2 ]} + 0(g«) 


-f^‘i*Ci(/S, + 4J 1 ,/ 2 ,)+0( 9 «) 

3o 4 

d A C 2 A (d - 1 )(d/ bb u + 4J la / b un + 2 J 3b ) + 0(/) 
3o 4 

-^d A CU-9(d - 1)/& - 12■/!»[/,»„„ + (2d - 3)/ 2 J 
+4(—2J 3b + 3J 3e — J 3 j + t/i a J 2a )} + 0(g 6 ) 


(A7) 

(AS) 

(A9) 

(A 10 ) 

(All) 


x{C A [(3d - 10)/ s b n - (lOd - 12)/ 2a - 8(d - 2)/ 2b - 2(d - 2)(d 2 -6 d + 6)M 2 ] 
+*S'f[— 4(d — 4)/g Un + 16(d — 3 )/ 2c + 32 / 2c j + 8(d 2 — 6d + 6 )/i& 2 ]} 


A^aC 2 


16 


Jl 


-2(d-3)^ + U la 


3d — 5 T d — 3 T „ t n 
-J 2 & + 4——-J 2c + 8 J 2( j 


d — 2 


d- 2 
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1 5d — 11 T 

J 3a + 2— - — J 3b + 


d — 2 
16 


d- 2 Jsg 


d — 2 
3d — 5 ' 

8 Y^2 J3h , 


d — 2 

+ 0(g 6 ) 


d3c — 64 J 3e + 64J 3 f 


-A 2 ‘F,„ = (C + 4J la / 2a ) + O(Z) 

-/i 2 e F n = 2d F "£ In P 2 

un 

72 „2 


-/i 2 ^ 0 = 


d A *S' F [3(d - 2)(/f - 2bifi) - 6(d - 2) J la A - 3M 2 (A‘ - 4A‘) 

O 

-(d-2)M 2 J la / 2 ] + 0(/) 


y e F p = ^d A (2S 2F - C' A S' F )(d - 2)[2(d - 4)/^ - (d - 6)/ b ff all + 8(6, - A)/^ 

+4Ji a / alm + 16J, a / 2a ] + 0(g 6 ) 

= 9 4 4 A S 2F {(d - 2) 2 [2ff 3 - - (6! - /,) 2 / 2 ] 

+2(4 — 2) J la [-/ S f un — 4I 2e + (4 - 4)(6, — /i)/ 2 ]} + 0( 9 6 ) 

-9 2< A = 9 4 4 A C A S F {(4 - 2)[/ b “ n + 2(6, - /,)/!] 

+4<7, a [(4 — 3)/2c + 2 / 2( i + (4 — 2)/ 2a + (4 — 3)/,6 2 ]} + 0( 9 fi ) 


(A 12) 

(A 13) 
(A 14) 


(A15) 


(A16) 


(A17) 


(A 18) 


The sum of those parts ill the contributions above that lead to (and potentially 9 s ln 9 ) 
terms as e —i 0 is 

- ^ 2e F (5) = -^d A M 2 J la {6(Z 2 - 1) + g 2 [(d 2 - 13d + 58 )C a 6 2 - 4(d - 2)P F / 2 ]} 

12 

,4 


2 y (d - 2)(d - 4)d A P 2F J la (6i - A)/ 2 

3d-5 


4 , „ 2 f d — 3 J ba 
+d d A C A A ---+ Aa 


8 M 2 
3d-5 
16(d — 2) 


4(d - 2) 


d — 3 

d 2 a H-;-yd 2c + 2 J 2 d 


d — 2 


7 (4 — 3) (3d — 5) 7 

Pa ~i-XT - ;-XT-Pb 


8 (d — 2) 


4 J 3e + 4 J 3f + 


-J- 


3d-5 


d — 2 3g 2(d — 2) 


2d-3 
2 (d — 2) 

d 3 h + d 3 i 


J- 


3c 


(A 19) 


where should be used up to order g 2 . Using the identities of Appendix [B] we can simplify 
this expression and get 


- h 2 ^ (5) = —-d A M 2 J la {6(Z 2 - 1) + g 2 [(d 2 - 13d + 58 )C A b 2 - 4(d - 2)P F / 2 ]} 
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. (A 20 ) 


-2 g\d - 2 ){d- 4)d A S 2F Jia( 6 i - /i )/ 2 


+g A d A C\ 


3d 2 -24d + 37 (c? — 3) s (3c? — 14) 

8 (c? — 7)(c? — 5 ) 2 3a + 4(d-5)(d-4)(2d-9) Jsb 
(d 5 - 23d 4 + 199c? 3 - 809c? 2 + 1548d - 1132)(d - 3) J 3 a ] 
8 (c? — 7)(c? — 5) 2 (c? — 4 ) 2 M 2 J 


Note that for this term the cancellation of overlapping double-frequency sum-integrals still 
holds outside of c? = 4. Using the results of Appendix [B] it is further easy to see how the 
1/e terms associated with the scales T and gT cancel separately so that for e —> 0 no g 5 In g 


term arises in F. 


APPENDIX B: BASIC INTEGRALS 

Here we give the definitions for the integrals appearing in our derivations. In the next 
subsection, we first provide the definitions of the sum-integrals which have been evaluated 
in Ref. jf|] and give the results for those that are relevant for the g 5 term of the free energy 
density. Then we define and give the results for five additional two-loop sum-integrals which 
appear in the result of individual diagrams but cancel each other after summing up the 
diagrams. In the second subsection, the definitions of and results for the three-dimensional 
integrals arising in the g 5 evaluation are given. They are evaluated in Appendix |7]. 


1. Some Sum-integrals 


Here is a list of integrals evaluated in Refs. (14]||. One-loop integrals are 
1 (2vr Tf~ 2n ( n 2 V T(n - § + e)C(2n - 3 + 2e) 


bn. — 


fn = 


p2r, 


87 T 5 / 2 \ nT 2 


T(n) 


f{p} P 2 
The relevant cases are 

T 2 


= (2 


2 n— 3+2e 


1 )b n . 


(Bl) 


bi = 


bo = 


1 + 2e (1 + In y 

12 V 4ttT C(—1), 

(- + 2 7E + 2hi-^) +0(e) 


+ 0(e 2 ) 


(4vr) 
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h 

/a 


T 2 
~ 24 
1 

(4^ 


1 + 2 e 1 - In 2 + In 


4 vrT 


+ 


C'(-i)' 


C(-1) 

(- + 27 b + 4 In 2 + 2 In + 0(e). 


+ 0 (e 2 ) 


Two-loop sum-integrals are 


^ “ Tpq P 2 Q 2 (P + Q ) 2 


/ f = ^ 
sun — 


l 


% {Q} P 2 Q 2 (P + Q ) 2 


= 0 
= 0 . 


(B2) 


(B3) 


Three-loop integrals are 

J bb = 'sf' _1_ 

bal1 jpqk P 2 Q 2 K 2 (P + Q + K) 2 

J bf = hf _1_ 

bal1 Tpq{k} P 2 Q 2 K 2 (P + Q + K ) 2 

/ba11 = f{P}{Q}{K } P 2 Q 2 K 2 (P + Q + K) 2 

H 9_ • K 

3 7|P}Qic P 2 Q 2 K 2 (P + Q) 2 (P + K) 2 ‘ 


(B4) 


Now we define some integrals that were computed in Ref. 

3po (1 — ^ go ) 


A\ = 


pq P 2 Q 2 (P + Q ) 2 


A h — ’sr Upo 
2 -TpqP* 


A[ = 


% 


Vo 


g 4 q 2 (p + q ) 2 


6 „ 


V 0 


p{q } P 2 Q 2 (P + Q ) 2 


At = sh f?E° 

2 - r P{Q} p 4 


% 2 


% 2 


g 4 g 2 (p + g) 2 


I (\C(\ — 


1 

P 4 


tr [n^(P) - n Att/ (o)] 2 , 


but not explicitly defined: 


(B5) 


where the values for the first four integrals may be found in the evaluations of hi, 82,63 in 
Ref. [Q and I qcd may be expressed in terms of / bb d , I q ^. d , and /q Cd there. 

The following five two-loop sum-integrals appear only in individual diagrams but not 
in the final result for the free energy density. They can be evaluated using the methods 
introduced in f|]. Here we only give the definitions and the values of these five integrals. 
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ha = 
hb = 
he = 
Pd = 

he = 


^ ^_ 

Tpq P 4 Q 2 (P+Q) 2 


g 0 2 


12(4vr ) 2 

T 2 


pq P 4 Q 2 (P+Q) 2 
Po 

mi p 4 g 2 (p+g) 2 

go 

mi p 4 g 2 (p+g) 2 

Pogo = _ 

{p}q p 4 g 2 (p + g) 2 ~ i6(4vr) 2 


+ 0(e) 


12(4tt ) 2 

T 2 

12(4tt ) 2 
T 2 


^ + I + I' 7b+ I 1112+3 In 47rr 


A 4 


3C' 

+ - 


- 1 ) 


3 7 3 5 3 

8 t + 8 + 4 7 ” + 4 ll2+ 2 11 4irT 


A 4 


2 C(-l) 

3 C ; (~l) 

4 C(-l) 


+ 0 ( 6 ). 


+ 0(e) 

+ 0 ( 6 ) 

(B 6 ) 


2. Three-dimensional Momentum Integrals 


Here is a list of our basic three-dimensional momentum integrals. We use dimensional 
regularization to control both the ultraviolet and the infrared divergences. Therefore, “three- 
dimensional momentum integrals” really means integrals in 3 — 2e dimensions. The steps 
for computing these integrals are provided in Appendix 0. 


•A i a. — 


J p 0 


p P 2 + M 2 
TM 


47T 


1 + 2 In 


jj_ 

2 M 


+ 1 e + 2 In 


o . 7T 

2 r + 4 In + 4 + — 


ha 

hb 

he 

+2d 

+3a 


ho ho 


2 M 

rp 2 


2 M 


+ 0 ( 6 3 ) 


1/2 1 

PQ (P 2 + M 2 )(g 2 + M 2 )(P + g) 2 (4vr) 2 L4e + 11 2M + 2. 


ho ho 


T 2 


1 , n l 

-h In-b In 2 H— 

4e 2 M 2 


8 (4n) 2 
1 T 2 


PQ (P 2 + M 2 )g 2 (P + g) 2 (4vr) 2 

M 2 h po h go _ 1 T 2 

p Q P 4 (g 2 + m 2 ) [(p + g) 2 + m 2 ] 

_ MH Po 5 qo _ = _ 

PQ (P 2 + M 2 ) 2 (g 2 + M 2 )(P + g) 2 4(4vr) 2 

_ ho ho h 0 _ 

pqk {P 2 +M 2 )(Q 2 +M 2 )(K 2 +M 2 )[{P+Q+K) 2 +M 2 } 
T 3 M 


+ 0 ( 6 ) 
+ 0 ( 6 ) 


+ 0 ( 6 ) 


+ 0(e) 


(4tt) 


- + 61n-^--41n2 + 
e 2 M 


+ 0(e) 


+3b = 


ho ho ho 


pqk (P 2 +M 2 )(g 2 +M 2 )/i 2 (p+g+/i ) 2 


T 3 M 

2(4tt ) 3 


61 ll 2 S + 


+ f 18 In 2 -4- + 48 In + 52 + 25?r 


2 M 


2 M 


12 


+ 0(e 2 ) 
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J 3c 
^3d 

hie 

<hf 

hs g 


M 2 6 po 6q 0 S ko 


T 3 M 


pqk K 2 (P 2 +M 2 )(Q 2 +M 2 )[(K+P) 2 +M 2 ] [(K+Q) 2 +M 2 } (4tt ) 3 

_ M2 W k0 _ = ™£ 21n2 + 0(e) 

pqk K 2 (P 2 +M 2 )(Q 2 +M 2 ){K+P) 2 (K+Q) 2 (4tt ) 3 

_ M 2 5 Po 5 qo 5 k „ __ T 3 Mtt 2 

pqk (Jl 2 +M 2 )(P 2 +M 2 )(g 2 +M 2 )(A>P) 2 (/l+g ) 2 (4tt) 3 12 


In 2 + 0 (e) 


+ ^( e ) 


^3h 

h i 


_ M\ 0 6 qo 6 ko _ = pgw; /V _ i\ 

pqk (K 2 +M 2 ) 2 (P 2 +M 2 )(Q 2 +M 2 )(K+P) 2 (K+Q) 2 (4tt ) 3 1^24 4 / U 

_ M\ 0 5 qo 5 ko _ 

pqk K*(P 2 +M 2 )(Q 2 +M 2 )[(K+P) 2 + M 2 ][(K+Q) 2 +M 2 ] 

T 3 M ( 1 In 2\ . 

Wl24 + ljj +0(£) 

M 2 S„S„S ta 


pqk K 2 (P 2 +M 2 )[(K+P) 2 +M 2 ]Q 2 (K+Q) 2 

M\ 0 8 q0 S k0 


T 3 M 

(47r) 3 


1 3, a 1 

- 1 — In- 

8 e 4 2M 4 


+ ^( e ) 


pqk (P 2 +M 2 ) l(P+K) 2 +M 2 ] (Q 2 +M 2 ) [( Q+K) 2 +M 2 ]K 2 (P+Q+K ) 2 

T 3 M /l In 2 \ . 

+ 0 (e) 


(47r) 3 \^4 4 


= 


p • g) 


pqk (. K 2 +M 2 )(P 2 +M 2 )(Q 2 +M 2 )(K+P) 2 (K+Q) 2 


T 3 M /vr 2 l\ _, , 

(top (12 + 1) + 0(e) 


(B7) 


APPENDIX C: EVALUATION OF BASIC INTEGRALS 

Here we will first evaluate J la , J 3a and J 3 b and then express all other three-dimensional 
integrals appearing in the diagrams (a)-(r) in terms of these three. 


1 . Jia 


J ia may be evaluated as 


•ha = T/J 


2e 


d 3 


—2e 


P 


= TM (—) 
(2n) 3 ~ 2e p 2 + M 2 KM) 

2 e roo J 1 M 


2 e 


™(s) 'C dss ^’ = 


d 3 2e p 
( 2 vr) 3 - 2e 
M 2 \ ~ e 


(47r) 3 / 2 \47r/x 2 


dse~ s ^ 2+1) 

Jo 

r(-| + e) 


TM 

47T 


1 + 2 (In — + l) e + (2 In 2 — + 4 In — + 4 + — ) e 2 

1 2 M J \ 2 M 2 M 4 J 


+ 0(e 3 ). (Cl) 
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2. J 3a and J 3b 


In dimensional regularization, we have 


J 3a — 2 J 3 b — X 5, 


if fco |?p V( p2 + ^ 2 )[(p + a') 2 + m 2 ] p 2 (p + /i) s 


I 2 


''PO 


^po 


(C2) 


We are now going to evaluate this difference using a method similar to the one applied in 
the appendix of the second reference of [ITU . 


Since J 3a — 2 J 3b is both infrared and ultraviolet finite, we can go to three-dimensional 
coordinate space to get 


J 3 a — 2 J 3b — T 3 / d 3 r 


,-2 Mr 


1 2 


(47rr) 2 (47rr)" 


+ 0(e) , 


where we have used the Fourier transform 

f d 3 p e l P' f 


J (27t) 3 p 2 + m 2 47t r 

for 77i = M and rn — 0. Integrating by parts and using the identity 


roc 1 h 

/ dr-(e~ ar — e~ br ) = In - 
lor a 


(C3) 


(C4) 


(C5) 


gives 


T 3 M 

^ 3 a — 2J 3 b = ——^-4 In 2 + 0(e). 


(4tt) 


J 3 b may be evaluated as 


J 3 b = T 3 M 


= T 3 M 


6e 


M7 J (2vr) 3 - 2e 

p 


d 3 2e p f d 3 2e q 


d 3 2e k 


(C6) 


(27r) 3_2e g 2 (g + p) 2 J (27r) 3_2e (fc 2 + 1) [(fc + p) 2 + 1] 


(i+ e ) 

r d 3 2e p r 

h-^ 

CO 

1 

to 

Os 

( 27 r) 3 — 2e Jo 


da[a(l—ct)p ] 


21-1/2-e, 


dp[p(l-(d)p 2 +l] 


—1/2—e 


= T 3 M f -^) 6e ^ - B ^ 2 6,2 f 1 dp [°° dp 2 (p 2 )~ 2e \P(P~P)p 2 +A i/2+e 

\M) (4vr) 3 - 2e (4tt)3/2-T (|— e ) Jo h Jo y ; ^ J 


= T : 


! M( — ) / 2 ' -^- 2 ; N P (l-2e,-^+3e) P (2e, 2e) 

\M) (47t) 3 - 2e ( 4 7 t) 3 / 2—e T (|—e) V 2 > y ’ 
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(47r) 9 / 2 \47T/i 2 

T 3 M Tl 

2(47r) 3 e 


r i 

{i+ e ) 

i r i 


- e ) 

i 2 r 

(-i+3e)r(2t) 2 

ri 

(i- 

-0 

r(4e) 


61 n 2 ^ + 8 ) + ( 181 n 2 2 ^ + 481 n 2 ^ + 52 + ^f) e 


+ 0 ( 6 2 ). 

(C7) 


For step two above, we have used the result 


d 3 


—2e, 


r(| + c) r 1 


da[a(l — a)p 2 + m 2 ] 


21 - 1 / 2-6 


(C8) 


J (27r) 3_2e (p 2 + m 2 ) [(p + g) 2 + m 2 ] (47r) 3_2e Jo 

obtained by Feynman parametrization, for m — 0,1. For step three, the identity 

f dxx a ~ l {l — x) 6_1 = B(a, b ), (C9) 

Jo 

where B(a,b) is the Beta function 

B(o ’ 6 ) = rrT$’ (C10) 

f [a + b) 

as well as the surface area of the d -dimensional unit sphere, 2n d / 2 /T(d/2), have been used. 
For step four, we used the identity 


fOO 1 

dx 77 ——r—r = B(a, b ) . 


Jo ^(l + x) a+b 
Combining the results (|C6|) and (|C7|) gives the value for J 3a . 


(Cll) 


3. Identities for Three-dimensional Integrals 

Here are some identities for (3 — 2e)-dimensional integrals. It is easy to derive them and 
we will give a sample proof in the following subsection. 

d-3 Jl 
2a 2{d — 4) M 2 

J 2c = —j - ~J2d 

d — 6 

J 2d = J 2a 
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d 3 c — 


4(d-5) 
d 3 d = —; -— d 3 b — 2Ji a J 2b 


[(3d — 11) J 3a — 4(d — 3) Ji a </ 2a . 


•1 3 e — 


d — 3 
1 


-[(3d — ll)</ 3b — 2 (d — 3) Ji a J 2a 


4(d — 4) 

d 3 f = _ gj [(^ — 13) J 3e — 2 (d — 3) Ji a J 2 d] 

d 3 g = [ (3d — 13) J 3c + 4(d — 3)Ji a d 2c ] 

3d -11 

J3h _ 4(2d-9f 3b 

d 3 i = 777 —— [(3d — 13)(2J 3e — J 3c ) + 4J 3f ] 


J* = 


4(d — 5) 

1 


d 3 b + J 3 d — 4J 3e — 2Jia(2J 2a — J 2b ) — 


du 

Ad 2 


(C12) 


Putting all of them together lets one express all three-dimensional integrals in terms of J la , 
d 2b , d 3a and d 3b : 


d 2 a — 


d~3 J? a 

' 2a 2(d - 4) M 2 

t d -3 djj 

2c 4(d - 6) Ad 2 


d 2 d — 

d 3c = 


d~3 d 2 a 
4 Ad 2 
1 

_ 4(d — 5) 


(3d — ll)d 3a + 


2(d — 3) 2 df a 


d — 4 M 2 


7 - 3rf - n / 97 7 

^3d ~ —i -— J 3b — z Jla<J2b 


d 3e ~ 


d — 3 
3d - 11 


d 3 b + 


(d-3) 2 df. 


la 


4(d-4) 4(d — 4) 2 Ad 2 


T _ (3d-13)(3d-11) 7 (d — 3) 2 (d — 5) Jf a 

3f / 7 /( \ / o 7 r\ \ 3b 


8(d — 4) (2d — 9) 


T _ (3d - 13)(3d- 11) j , (d-3) 2 

^3g — ; "3a + 


8(d — 4) 2 M 2 
3d - 13 


16(d — 5)(d — 7) 
3d-11 

J3h “ ~4(2d-9) Jsb 

(3d - 13)(3d- 11) 
16(d — 5) 


8(d — 7) 


+ 


(d — 4)(d — 5) d — 6 


J, 3 


la 


M 2 


d 3i = 


1 7 4 7 ' 

d 3 a + 77- 7 - 7 3b 


d — 5 


T _ d 2 — lOd + 23 T 

3j A ( _1 A \ 1 _/ O \ 3 ^ 


2d-9 

7 3 

J la 


+ 


(d — 3) 2 (5d — 23) df. 


la 


8(d-4)(d-5) 2 M 2 


4(d — 4)(d — 3) 0D 4(d — 4) 2 Ad 2 


(C13) 
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Since in the diagrams, J 2 b and J 3d only appear through J 3] and therefore only in the com¬ 
bination J 3d + 2 J la J 2 b = (3 d — 11) J 3 b/(d — 3) [see Appendix [A] and Eq. ( |C12|) 1, we have not 
bothered to write down the evaluation of Jo b, although it is easy and can be done in general 
dimension. 


4. Proof of Identities for J 3c and 


As an example, here is the proof of the identity for J 3c in (|C12|) . The proofs of all the 
other identities proceed along the same lines with the exception of that for J 3 j, which is 
presented below. 


Using the shorthand 


T p 


2e 


d d l k 


(2jt) 


d-1 


(C14) 


we can write 




M 2 


kJpJq k 2 (p 2 +M 2 )(q 2 +A4 2 )[(k + p) 2 +Ad 2 }[(k + q) 2 +M 2 } 


M 2 


k, 


d 


M 2 


d-1 JkJpJq dki k 2 (p 2 +M 2 )(q 2 +M 2 ) [{k + p) 2 +M 2 ][(k + q) 2 +M 2 ] 


M 2 


d — 1 JkJpJq 
1 
d 


2 + 


4 k ■ (k + p) 


(.k + p) 2 +M 2 

—r i 4 J 3c + 2M 2 f j f 
— II JkJpJq 


k 2 (p 2 +M 2 ) ( q 2 +M 2 ) [(k + p) 2 +M 2 } [(k + q) 2 +M 2 } 

1 


(p 2 +M 2 )(q 2 +M 2 ) [{k + p) 2 +M 2 } 2 [(k + q) 2 +M 2 } 

1 


k 2 (q 2 +M 2 )[(k + p) 2 +M 2 ] 2 [(k + q) 2 +M 2 } 


1 . ,o dJ 3a , , ^2 d Jla 


d -5 
2 

d — 5 


--M 2 —— 
4 dM 2 


+ M* 


dM 2 


J‘ 


2a 


3d — 11 d- 3 \ 

"d 3a H -- - Jl&J2a, J 


8 


(CIS) 


where for the second equality we have used the identity (dki/dki) = d — 1 and integrated 
by parts and where for the last step it has been used that by dimensional considerations 
J 3a oc M 3d ~ u and Ji a oc M d ~ 3 . 


Finally, here is the proof of the identity for J 3 j. Noting that 


























p 


Ip ( p 2 +M 2 )(k + p ) 2 


(C16) 


is parallel to k, we have 


■Jm = 


p-q 


IkJpJq (k 2 +M 2 )(p 2 +M 2 )(q 2 +M 2 )(k + p) 2 (k + q) 2 

f f f _ (k ■ p)[k • q) _ 

kJpJg k 2 (k 2 +M 2 )(p 2 +M 2 ){q 2 +M 2 )(k + p) 2 (k + q) 2 


j 3 

^3b + J'id ~ 4J3e ~ 2Ji a (2J 2 a “ ^2b) ~ 


(C17) 
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